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Chapter 1 
Introduction 



One of the most important and unresolved questions of canonical gravity is the so called 
"problem of time" . In quantum theory time is not regarded as an observable in the usual 
sense, since it is not represented by an operator. Time is rather treated as a parameter 
which, as in classical mechanics, describes the evolution of a system. In the time-dependent 
Schrodinger equation, for example, time occurs only as parameter and not as operator. In 
the conventional Copenhagen interpretation the notion of a measurement made at a par- 
ticular time is a fundamental ingredient. An observable is something whose value can be 
measured at a fixed time. 

In general relativity the role of time is very different. The four dimensional spacetime, 
i.e. a four dimensional manifold with a Lorentzian metric and an appropriate topology, can 
be foliated in many different ways as a one-parameter family of spacelike hypersurfaces. 
Each of these parameters can be viewed as a possible definition of time. There are many 
such foliations and there is no way to select a particular one to consider it as 'natural'. 
Hence, in general relativity, time is not a fixed quantity, as it is in quantum theory. Time 
is dependent on the chosen foliation of the manifold and is closely linked to a fixed choice 
of the metric. Such a definition of time does not give any hint how it should be measured, 
and is so a non-physical quantity. 

The gravity group of the institute of theoretical physics deals with asymptotically flat 
gravitational models. This kind of systems seems to remove the "problem of time" in the 
asymptotic, since the parameter of time can be chosen in a unique manner. The dynamics 
of the system can be determined, and one can attempt to study theories like the scattering 
one in a quantum theoretical way. The hope is to understand the structure of the theory of 
such models in order to be able to state something in the direction of a quantum gravity. 

In this diploma thesis such an asymptotically flat model is investigated. The model I 
have chosen in this paper is pure gravity - no matter fields are present - with cylindrical 
symmetry. The symmetries of the cylindrical gravitational wave help to solve the Ein- 
stein equations and thus, to know the metric which describes spacetime of the treated toy 
model. The cylindrically symmetric gravitational wave, also called Einstein-Rosen wave, 
was discovered by A. Einstein and N. Rosen [1]. To give an idea what gravitational waves 
are we cite [2]: Just as one identifies as water waves small ripples rolling across the ocean, 
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so one gives the name gravitational waves to small ripples rolling across spacetime. In 
the seventies of the last century this field theory was again picked up by K. Kuchar [3]. 
The model was solved in the canonical ADM formalism and afterwards quantized. Twenty 
years later, C.G. Torre [4] found a complete set of observables of the cylindrically symmet- 
ric gravitational wave. 

This work is mainly based on the contributions mentioned above. In chapter two the 
general cylindrical symmetric metric is studied and mathematical simplifications of it are 
discussed in detail. Then the Einstein-Rosen metric is introduced and the corresponding 
Einstein equations are derived. In the third chapter, we determine the solution of the 
Einstein equations and check the boundary conditions the solving function is supposed 
to meet. In the following chapter the ADM formalism is presented, firstly for the gen- 
eral spacetime then for the cylindrically symmetric one and finally for the Einstein-Rosen 
wave, which is the interesting case for our consideration. The fifth chapter is devoted to 
the derivation of the dynamical variables by solving the Hamilton equation. Moreover it 
is shown how in Torre's work the set of observables for the cylindrical gravitational waves 
is obtained. In the last chapter the condition is investigated, which phase space functions 
have to fulfill in order to be declared as observables. Then the Poisson algebra of the 
observables is studied and quantum theoretical conclusions to these quantities are stated. 
In the appendix a detailed derivation of the Hankel transformation is noted, as well as the 
orthonormality relation of the Bessel functions of the first kind. This relation plays an 
important role in the mathematical part of this diploma thesis. Lengthy computations and 
verifications have also been displaced into the appendix in order to preserve the overlook 
while reading. 

Conventions and notations: 

Spacetime has the signature S = 2, so the signs of the diagonal elements of the Lorentzian 
metric are (-, +,+,+). We use natural units that means K = c = IQ-kG =1. In order to 
distinguish clearly the coordinate R from the Ricci scalar R, we write throughout the coor- 
dinate in the italic style and the curvature scalar in the Roman style. The Greek indices run 
over 0,1,2,3 while the Latin indices only over 1,2,3, corresponding to the three space coor- 
dinates. Derivatives in general are denoted as follows: ^ f(T, ^) = > ~§Tvt(T, R) = f, TT . 
In particular the derivative with respect to the time t is denoted with a dot, the one with 
respect to the radius r with a prime. All these conventions hold unless otherwise specified. 
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Structure of the Cylindrical 
Spacetime 

2.1 Line Element with Cylindrical Symmetry 

A method to quantize a gravitational system is to foliate the four dimensional spacetime 
in a (3+l)-dimensional general relativity and subsequently to reduce it to an equivalent 
(l+l)-dimensional model, which is coupled to a massless scalar field. An example therefore 
is the cylindrically symmetric Einstein-Rosen wave. In this section we are going to study 
the line element for the mentioned kind of gravitational waves. 

We start with considering the line element, which defines cylindrical spacetime, in its 
most general form. Afterwards we will deduce from the general form, the line element 
in the Einstein-Rosen coordinates through conformal deformations and coordinate trans- 
formations. First of all, we take for convenience the cylindrical coordinate system to 
describe points in spacetime, namely the set of variables (t, r, tp, z) with t G (—00, 00), r G 
[0,oo), ip G [0,27r) and z G (—00,00). The most general line element with cylindrical 
symmetry reads: 

ds 2 = g a /3 x a x 13 = g o dt 2 + 2 g i dt dr + g n dr 2 + g 22 dtp 2 + g 33 dz 2 , (2.1) 
with obviously 

a; = t, x 1 = r, x 2 = ip, x 3 = z. (2.2) 

To get a cylindrically symmetric line element the coefficients of the metric g a p has to 
depend only on the coordinates t and r. As we are going to see, only in that case the line 
element is invariant with respect to changes of the coordinates ip and z. So the metric is 
a function of the two coordinates t and r, and of the form 



( 9oo(t,r) g 01 (t,r) N 

9w(t,r) gu(t,r) 

g 22 (t,r) 

\ g 33 (t,r) I 



(2.3) 
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where, due to the symmetry property of both indices of the metric tensor g a p, goi(t,r) = 
9w{t-i r )- After these settings, we now derive the symmetries of spacetime by analyzing the 
line element ( 12. ip . For this inspection we define the z-axis as the symmetry axis for the 
following four symmetries: 

1. Rotations around the symmetry axis z: ip — > ip' = ip + Atp. 

2. Translations in direction of the z-axis: z —>■ z' = z + Az. 

3. Reflections in all surfaces containing the symmetry axis: <p — > ip' — —ip. 

4. Reflections in all surfaces perpendicular to the z-axis: z — > z' — —z. 

If we substitute in the line element z with z' and ip with ip' the line element will not change. 
This is the reason why we have pretended g a p to depend only on t and r, otherwise the line 
element would change because of the substitutions. Seeing that, ds 2 is invariant under the 
four transformations noted above, just the four symmetries. The line element is symmetric 
to the following two dimensional transformation group: 

ip = ip + Aip (rotation around z-axis) (2-4) 

z — > z = z + Az (translation in direction of the z-axis) (2-5) 

All transformations of such a group are called isometries. Each of these isometries is 
generated by the related Killing vector field, for the cylindrical coordinates (£, r, <p>, z) the 
two orthogonal vectors 



= (0,0,0,1), (2.6) 



eg,) = (0,0,1,0) 

tp 
•to 

which are determined by the Killing equation 

g a /3,u C + 9vp C,« + 9au Q = (2.7) 

or equivalently 

d K g pa = 0. (2.8) 

The two isometries and consequently both Killing fields implicate conservation laws of 
the dynamic. Indeed the cylindrical gravitational fields are characterized by the existence 
of a two-parameter Abelian group of motions 4 G*2 with both, to each other orthogonal, 
hyper-surface-orthogonal spacelike Killing vectors 

2.2 Simplification of the Cylindrically Symmetric 
Metric 



In order to obtain solutions of the Einstein equations written for the general cylindri- 
cally symmetric metric ( 12. 3p . we have to simplify it without loosing physically relevant 
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properties. The method, which in our case is going to be applied, consists in deforming 
conformally a part of the four dimensional cylindrical metric ( 12.31) . and subsequently to 
find a coordinate transformation, which leads into an inertial system IS. Also we will see 
that the mentioned coordinate transformation brings the metric, dependent on the new 
coordinates, into diagonal form. 



2.2.1 Conformal Deformation 

In the beginning, we would like to start with defining the conformal deformation in space 
time as follows: 

g a p = IT 2 g' af3 . (2.9) 

The function Q = Q(t, r) is a real and positive scalar function. We consider now a conformal 
deformation in two dimensions. In this case the metric g\ K , with A, k — 0, 1 labeling (t, r), 
corresponds to a regular 2x2 matrix. From this it follows that 

i 

<?a k <? AK = ^ = $> A A = 2, (2.10) 
o 

where in section 2.2 A, k, p, a, /i and v run over and 1. Then we are going to show 
that there is always a conform deforming factor such that the metric g\ K can be expressed 
through a flat metric g' Xti . 



Theorem: Every two dimensional metric g\ K is conformally flat. 



Proof: We begin with affiliating the transformation law between the Ricci scalar of a 
conformally deformed metric, R[<7 A J, and the curvature scalar of the original one, R[(?ak]- 
From the definition of the conformal deformation (12.91) it follows: 



'Ak 



2 „Ak 



The Ricci scalar R is defined by 



R = g XK R AK = g AK R 



Ak dP 



Apft ' 



(2.11) 



(2.12) 



where R\ K is the Ricci tensor and R^„^ the contracted Riemann tensor. The transformed 



curvature scalar reads: 



XpK 



R' = g'*» R Att = g'^ R^. 



(2.13) 

If we put in (12.131) instead of g\ K consequently Q 2 g XK and write down the transformed Ricci 
scalar R', we get 



R =Vl 2 g XK 
Evidently, the Christoffel symbol 



8 F p 

U P L Ak 



F) v p -l r p r v —vpv u 



(2.14) 



rP A« = \ [duK, X + 9aX, K ~ 9Xk, u ] 



(2.15) 
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has to be transformed. Its transformation law is given by utilizing (12.111) . as done for the 
curvature scalar (12.1 2p . 

r'^ = \Vi 2 g p ° [(tt- 2 g aK ) ,a+ (^Va) , K -(n- 2 g XK ) „] . (2.16) 
If one performs the partial derivatives and the relation 

lnfi,„=(lnn),„=^ (2.17) 
is employed, then the transformed Christoffel symbol is obtained: 

r' a. = r p AK - in n, x 5 P - in n, K s p + in n, a g XK g p °. (2.18) 

This intermediate result is introduced in (12. 14j) and the transformation law for the Ricci 
scalar is given: 

RVaJ = n 2 R[g XK ] + n 2 (2 d K d K In Q + 2 T Kp In Q, p ) . (2.19) 

To this end we have used the fact that the covariant partial derivative of a scalar field 
(scalar function) is a covariant vector field. 

V p = $, p =d p $ = \nn, p . (2.20) 

From this it follows the repeatedly applied equation 

g XK d K V x = g XK (d K V x ) 

= g XK d K (V p g pX ) 

= 5;d K V p + V p g XK g pX , K 

= d K V K + V p g XK g pX}K . (2.21) 

Further more the subsequent expression is relevant also: 

2V p T x pX = g XK g XK , p V p . (2.22) 

In case that Q is a scalar field, the right hand side of the equation (12.191) can be noted 
more compactly. So, we apply the (two dimensional) Laplacian operator A onto a scalar 
field in the two dimensional spacetime. 

A0 = V k (s Ak V a 0) 

= V K V K 

= d K V K + T% K V p 

= d K d K $ + T K pK d p $, (2.23) 

whereby we know that the covariant derivative of a scalar field is equal to the covariant 
partial derivative of the same. From (I2.2UI) we see that $ = lnf2, and by putting it into 
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f)2.19p . we achieve that the transformation law for the Ricci scalar reduces to the following 
compact formula: 

R'[g' XK ] = Q 2 R[g Xfi ] + 2 ft 2 A In ft. (2.24) 

Now we set the two dimensional conformally deformed metric g XK to be flat. In this case 
the transformed Ricci scalar R'fg^J vanishes. From the transformation law ("12.241) it follows 
then: 

A In ft + \R[g XK ] = (2.25) 

This equation is the known generalized inhomogeneous Klein-Gordon differential equation 
for the massless scalar field $ = In ft. This kind of differential equation can be worked out 
by performing a change of variables to the double null coordinates: 

u = t + r, v = t-r. (2.26) 

In [3] it is shown that by the substitution to the double null coordinates the differential 
equation (12.251) possesses for all regular ft > a solution. So it exists throughout a scalar 
function ft such that the Klein-Gordon equation is fulfilled. Thus, there is a conformal 
deforming scalar function ft in such a way that the metric g\ K is expressed by a flat metric 

2.2.2 Coordinate Transformation into the Inertial System 

In the previous paragraph it has been showed that the (t, r)-part of the cylindrically sym- 
metric line element is conformally flat. By virtue of this property it can be expressed by a 
two dimensional flat metric. 

g XK dx x dx K = ft" 2 g' XK dx x dx K , (2.27) 

9oo = tt~ 2 goo, g' m = ft~ 2 #oi, 9n = &~ 2 g n . (2.28) 

The scalar function is chosen such that R'[g'\ K ] = holds, in order the metric g' XK to be 
flat. But by then, the four dimensional metric (12.31) has not been simplified really, since 
the number of its non- vanishing components did not decrease. Such a simplification would 
be an advantage to find solutions to the vacuum Einstein equations, as the number of 
vanishing components would rise. At this place it is again important to remind that the 
physical model has not to be trivialized by the possible elimination of some components 
of the Einstein equations. The issue is only to obviate the mathematical impediments by 
a suitable choice of coordinates. 

The theorem about the conformal flatness of two dimensional metrics just allow us 
to find a global coordinate transformation to an inertial system, as the conformally flat 
metric can be described by a flat metric g' Xft . Now the importance of the proved theorem 
appears evidently on the score of the suited simplification of the Einstein equations by 
reducing the number of non- vanishing components of the four dimensional metric (12. 3p . In 
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9ap(t,r) 



(2.31) 



the coordinates t and r of the inertial system the conformally deformed two dimensional 
(t, r)-part of (I2.3p . namely g Xlt , just transmute to the Minkowskian standard metric: 

i / . \ I ttt TTT i coordinate transformation w ,— _% / — 1 \ 

0ak(*,O = ( m m. J ^ »7ak(*,0 = ( Q 1 \ (2.29) 

Through the conformal deformation Q and the coordinate transformation (p we get from 
the conformally flat metric g\ K {t,r) a diagonal matrix, namely 

9Xn{%r) ^ > gx K (t,f)=n 2 (t,f) VxK (i,f) = (-Q 2 (2.30) 

The conformal deformation and the coordinate transformation simplify the exclusively on 
t and r depending four dimensional metric (12.31) significantly, as it is brought into diagonal 
form without changing anything on the physical aspect. The four dimensional metric reads 
newly: 

/ -Ci 2 o o o \ 

tt 2 
g 22 (i,f) 
V g 33 (i,r) J 

with the knowledge that #22 (i, f ) and #33 (?, f) are invariant with respect to the coordinate 
transformation. This fact is shown to hold in equation (I2.4ip if we are allowed to anticipate 
#22 (^ r) = R 2 e~^ and ^ 33 (i, f) = e^ in case of the Einstein- Rosen spacetime. 

At last we would like to add an illustration for a better intuitively understanding of 
the coordinate transformation into the inertial system. The starting point is a flat geom- 
etry, which is represented by the so-called flat metric g' Xn in (12.291) . However we probably 
would not assert on the basis of the treated metric that a two dimensional flat geometry 
is described, since for this case we rather would expect the Minkowskian standard metric. 
Nevertheless, we get zero by calculating the Ricci scalar of g' XK ; then we say the metric 
is flat. How should we understand intuitively this apparent paradox? To this end, we 
consider a flat disk at rest in an inertial system. With respect to the coordinates (£, f) of 
an inertial system, an observer A measures the length C = 2n for the circumference of the 
disk with the unit radius. Now we make the assumption that an other observer B sizes 
the disk with respect to a coordinate system rotating around the symmetry axis, which 
intersects the middle point of the disk. We emphasize that the rotating coordinate system 
is accelerated. The observer B also measures the radius to be unit, but a circumference 
C' smaller than 2n. This is an effect due to the Lorentz contraction. So the observer B, 
confident in his ratios, would not assert the measured surface to be a disk, but rather a 
paraboloid. Thus, he would confirm a curved geometry and consequently a two dimen- 
sional metric, which is different from the Minkowskian standard metric, presumably with 
four non-vanishing components. However the Ricci scalar of his metric is zero implicating 
a flat geometry. 

In a certain sense, he has applied for his gaging the "wrong" coordinate system, namely an 
accelerated one, which intuitively speaking, distorts the geometry. So it is recommendable 
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to undertake a coordinate transformation to the inertial system to see the geometry the 
"right" way. That the Ricci scalar of a flat metric vanishes with respect to all coordinate 
systems is not surprisingly, since a scalar is independent on coordinates. That is the reason 
why also observer B has to get a vanishing Ricci scalar for his metric. 

By means of this Gedankenexperiment, we see that the physical model is not influ- 
enced at all by a coordinate transformation. In fact the properties of the two dimensional 
geometry are not changed, it remains flat with respect to all coordinate systems. 



2.3 Einstein-Rosen Metric 



After the introductory sections, we shall consider more concretely the real subject of our 
study, namely the Einstein-Rosen gravitational waves. The utilized metric for this model 
is the Einstein-Rosen metric, namely 



( 


_ e r-^ 








\ 























R 2 e"^ 





\ 











e^ / 



^2.32) 



The functions V and ift are dependent on both Einstein-Rosen coordinates T and R ex- 
clusively. The complete coordinate system of course is eked by the other two spatial 
coordinates tp and z. As we notice, the Einstein-Rosen metric is cylindrically symmetric, 
since it is a special case of the general cylindrically symmetric metric (I2.3p . The metric 
( I2.32j) depends on the time coordinate Tg [0, oo) and on the spatial coordinate R G [0, oo) 
only. Just this requirement was the second requisite in section 2.1 a metric is supposed to 
possess in order to be cylindrically symmetric. Thus, the Einstein-Rosen metric inherits 
all properties of cylindrical metrics and which are listed in section 1.1. 

In the following it is showed, how one gets from the general cylindrically symmetric 
metric (12.31) the Einstein-Rosen metric, which is the important special case for our purpose. 
If we set in the general metric ( 12.31) 



g 22 (t,r) = R 2 e-^ r)=e* 
Hwe remark that only the two dimensional (t, r)-part 



9oo{t,r) g i{t,r) 
9oi(t,r) gn(t,r) 



g\ K (t,r) = 

has to be brought into the diagonal form 



(2.33) 



(2.34) 



9x K R (T, R) 



-e r "^ 







S—l> 



-1 
1 ' ' 



(2.35) 



1 Compare with equation (2) in section II of [3] 



10 



Chapter 2. Structure of the Cylindrical Spacetime 



where A and k run again over and 1. In the original coordinates (t, r, tp, z) the functions 
R and ip are dependent on t and r. In Kuchar's paper [3] both functions are defined by 
the Killing vectors ( 12. 6ft : 



ij(t,r)=\n(^ z) Z a{z) ). 



(2.36) 
(2.37) 



If the right hand side of both equations ( 12.361) and ( 12. 37ft are calculated by means of both 
Killing vectors and the metric 



f 


goo 


001 








\ 




0oi 


0n 


















R 2 e"* 







\ 













/ 



(2.38) 



then the functions are got. Now we bring g\ K (t, r) (12.341) into diagonal form. The metric 
is two dimensional and dependent only on t and r. After we have learned in section 2.2 
how this kind of metric can be simplified, the following plan appears to be suited: 

1. We find a scalar function Q such that the metric g\ K can be expressed by a flat metric 
g' XK . To this end the Klein-Gordon differential equation for the massless scalar field 
Q (ESD must be solved. 

2. By a general coordinate transformation (p, a coordinate system (£, f) of an inertial 
system IS is found, wherein the flat metric g'\ K (i,f) assumes the Minkowskian stan- 
dard form. 

The result of this procedure is the metric which is conformally flat in the (t, r)-part. This 
metric is noted in the second section of Kuchar's paper [3] , 



/ 













\ 

























R 2 e-^ 







V 













/ 



(2.39) 



for the case that the solution of the Klein-Gordon differential equation f!2.25|) is the subse- 
quently conform deforming factor: 



Q 2 



[2A0) 



7 = 7(t,r), ifj = ij(i,f) = ij(t,r). (2.41) 

The functions R and i\) do not change because of the coordinate transformation ip, since 
they behave like scalar^, as one can take from the equations (1 2 . 3 6 f) and (12.371) . 

2 Now the provision we made for the metric (|2.31j) is justifiable in virtue of the behavior of R and if) 
with respect to coordinate transformations. 
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We shortly focus on the Einstein-Rosen coordinate R. To develop a more concrete 
vision, as also to convince us that the Einstein- Rosen coordinate is a kind of radius in the 
cylindrical space time, we add a short illustration. We concentrate on a two-dimensional 
cylindrical surface i=const., r=const. around the axis of symmetry. Let be the height of 
the cylindrical surface Az = 1. Thus the area of the described surface is 2irR, the same 
as the area for the surface of a cylinder with height one and radius R in Euclidean space. 
So we can identify R in (I2.3ip with a kind of radius in the cylindrical spacetime. To make 
sure that this quantity is spacelike, a special regard should be given to the vector R, u (£, r), 
which can be spacelike (R, U R, U > 0), timelike (R, u R, u < 0) or lightlike (R, u R, u = 0). In 
order to deal with the model of cylindrically symmetric gravitational waves, we are forced to 
choose the vector R, u (t, r) spacelike everywhere in spacetime. Of course, it would be quite 
senseless to speak of a timelike or a lightlike radius vector, at least if the coordinate R is 
supposed to form with the timelike Einstein- Rosen coordinate T an orthogonal coordinate 
system. 

2.3.1 Conformal Coordinate Transformation 

Comparing the metric g a j3{i, r) with the Einstein-Rosen metric g^(T, R), we notice that a 
further simplification of the first metric is possible. By a coordinate transformation $ into 
the coordinate system (T, R, tp, z), the function R(i, f) can be converted into a coordinate, 
namely one of the Einstein-Rosen coordinates. Looking at both metrics above, we remark 
that besides a further conformal deformation in the (t, r)-part 

9x* = 9t R , (2.42) 

just a change of coordinates transports the one metric in the other (please notice A, k = 
0, 1). Then the functions T(t,r) and R(i,r) adopt the role of coordinates T and R in the 
orthogonal Einstein- Rosen coordinate system (T,R,ip,z). While in these coordinates the 
function ip(t,r) will not change as it is a scalar, the function ^(i, f), on the contrary, will 
change due to the coordinate transformation, according to the relation, 

j(t, f) -> r(T, R) = 7 - In (R' 2 - T 2 ) , (2.43) 

where the primes are in this case derivatives with respect to f. Since the change of coordi- 
nates influences the time and the radius, we require the following connection between the 
(t, r)-parts of the considered metrics: 

- dt 2 + df 2 = Cl- 2 (-dT 2 + dR 2 ) . (2.44) 

Thereby it is not about a general coordinate transformation, but about a so-called confor- 
mal coordinate transformation. From this requirement we are going to get an additional 
condition on the coordinates T and R. We compute the total differential of the functions 
T(i, f) and R(i, f) and set them in the right hand side of equation (I2.44p after multiplying 
with fl 2 . So we get: 

n 2 (-dt 2 + df 2 ) = - (Tj -R 2 t )di+ (R 2 f -T 2 )df 2 

- T d T, f dtdf + R d R, r dtdf. (2.45) 
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The equation (12.451) holds only if one of the following two settings is applied: 

T,t = R,f , T, r - = R,t (2.46) 

T,t = —R,f i T, f = —R,t- (2.47) 

We have assumed that the order of the partial derivatives is commutable, what for coor- 
dinate transformations surely apply. Derivating partially the first equation in (I2.46P by i 
and the second in (12.461) by f, and then subtracting the second result from the first, we get 

□ T(t, f) = T, E -T, F r = 0. (2.48) 

Derivating partially the first equation in (12.461) by f and the second in (I2.46P by i, and 
then subtracting the first result from the second, we get 

□ R(t, f) = R, rt -R, ff = 0. (2.49) 

So, the Einstein-Rosen coordinates T and R not only build by construction an orthogonal 
system and are timelike respectively spacelike, but also they have to be harmonically 
with respect to the coordinates i and f. Due to the invariance of R regarding coordinate 
transformations, they are harmonically with respect to the original coordinates t and r, 
too. Of course the coordinate T is also a scalar and therefore invariant to an arbitrary 
change of coordinates, as the quantity follows by integrating the equations ( 12.461) or (12.471) . 

In order to verify whether the function R(i, f) is harmonic, we write the Einstein 
equations for the metric (I2.39j) : 

Ga/3 = — \ = (2.50) 

By performing the Ricci tensor K a p and the Ricci scalajf) R for the metric g a p, we obtain 
the following, for our consideration relevant, components: 

G 22 = R(l,tt -7,fr ) + \ R(lp,l ~1p, 2 r ) = (2.51) 

G 33 = R{lM -7,rr ) + \ R{H -t/v- ) + 2(R,u -R,ff ) = (2.52) 

Equation (I2.49j) just appears by subtracting G 33 from G 2 2, 

R,tt~R,fr — 0. (2.53) 

We assume the derivative of the Einstein-Rosen coordinate R to be by construction a 
spacelike vector (i?, y R, u > 0) and the coordinate itself to be perpendicularly on the time 
coordinate T. Then R is fit to be Einstein coordinate, since it is showed it to be harmonic 
with respect as well to the variables (i, f) as to {t,r). At last we would like to note both 
Einstein- Rosen coordinates T and R. To this end we solve both D'Alembert equations 
( I2.48P and (I2.49P for one time dimension and one space dimension. This sort of differential 

3 In the subsequent section the way how the tensor and the scalar are performed is explained in a more 
detailed way. 
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equation can be worked out by writing them in double null coordinates, i.e. the following 
substitution of variables is suggested: 

u = i+f (2.54) 
v = i-f. (2.55) 

If we note the differential operators in terms of the new coordinates, then we get for the 
time and radius derivatives 

d du d dv d d d . . 

di di du di dv du dv ' 

d du d ^ dv d d d (2 57) 

df df du df dv du dv 7 

and putting them in (12.481) and (I2.49p . they yield the D'Alembert differential equations in 
the double null coordinates: 



d 



2 



-T(u,v) = 0, (2.58) 



dudv 

^ = 0. (2.59) 

So, the general solution of the differential equation (12. 59f) can be found. For the derivative 
with respect to v the following is compelled to be set: 

^R(u,v)=f(v). (2.60) 

Only in case that the derivative with respect to v yields a function depending on v exclu- 
sively, the following derivative with respect to u vanishes. In that case the D'Alembert 
equation (12.591) is fulfilled. The indefinite integration over v of equation (12.601) yields the 
general form of the solution for the incoordinate: 

_ d _ _ f 
dv—R(u,v)= / dviiv) 
dv J 

-> R(u,v) = F(v) + G(u). (2.61) 

Analogously the general solution for the time coordinate T reads: 

_d _ _ f 
dv—R(u,v)= / dvh{v) 

T(u,v) = E(v) +K(u). (2.62) 
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As we know, both Einstein-Rosen coordinates are linked to each other by the relations 
( I2.46P or (12.471) . From the first identity in (12.461) the connection in double null variables 
follows: 

m <9G dF dK dE . 

T,;= R,, — > j- + Tp = J- + J- ■ 2.63 

du dv du dv 

From the second identity in (12.461) we get a similar expression performing in the same way 

_ dG dF dK dE , . 

T, f = R,t — > _- — = _ + —. 2.64 
ou ov ou ov 

Adding both identities written in the double null coordinates we obtain 

dG dK 

du du ' 

-> G(u) = K(u), (2.65) 
and by putting this result in ( 12. 63ft or ( 12.641) trivially 

dF dR 

dv dv 

-> F(v) = -H(«). (2.66) 

We remind on the fact that in (I2.65P and (I2.66P functions depending on u respectively v 
are not allowed to occur due to the indefinite integrations. Otherwise the condition 

JLJLr(u,v) = f(«) = 0, (2.67) 
du dv du 

i(v) = -jj=F(v), (2.68) 

would be violated. Replacing the double null variables by the original coordinates i and f 
the Einstein-Rosen coordinates read: 



T(t,f) = G(t + f) -F(t-f), (2.69) 
R(i,f) = G(i + f)+F(i-r). (2.70) 

As we have seen, the conditions ( 12.461) and ( 12.471) . which arise from the crucial requirement 
the transformation between (t, f) and (T, R) to be conformal, is authoritative. If we con- 
struct the coordinate R(i, f) then the coordinate T(i, f) emerge directly from the relations 
( I2.46P and ( 12. 47ft . It has been showed that the function R(i,r) is actually harmonic by 
solving the Einstein equations written for the metric (12.311) . Rather due to the relations 
(12.461) and (12.471) it is also proved that the time coordinate T(t, f), which corresponds to the 
Einstein- Rosen coordinate R(t,f), is also harmonic. Therewith the differential equation 
(g3BD is fulfilled, too. 



2.4 Einstein Equations in the Einstein-Rosen 
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2.4 Einstein Equations in the Einstein-Rosen 
Coordinates 

The model of the cylindrically symmetric gravitational waves does not include the presence 
of matter fields. Thus, we consider a general relativity without matter fields, that means 
a vacuum spacetime. So, we are interested in the vacuum Einstein equations 

G a /3 = R-a/3 — \ R ga/3 = 0. (2-71) 

In this case the energy-momentum tensor T a p is zero and the Einstein equations 
G a /3 — KT a /3 = reduce to the vacuum equations (I2.7ip . 

We now would like to calculate the vacuum equations (I2.7ip in the Einstein-Rosen 
coordinates system, since, as we are going to see, they assume a very convenient form 
therein. So, the metric g a p is obviously the Einstein- Rosen metric (12. 32ft . In virtue of 
the cylindrical symmetry, only six components of the Einstein tensor (I2.7ip do not vanish, 
namely: 

Goo = R-oo ~~ | R #00 Goi = Roi (2.72) 

Gio = Rio Gn = Rn — | R #11 (2.73) 

G 22 = R22 — \ R #22 G 33 = R 33 — \ R # 33 . (2.74) 

Because of the symmetry of the Ricci tensor R Qj g in its indices, the components Goi and 
Gio are identical, and so one of them can be neglected, for instant Gio- Now only five 
components are relevant for the further computation of the Einstein field equations. In the 
apposite appendix C the various non vanishing components of the Ricci tensor 

iW = = d » - dp + 1%. - ( 2 - 75 ) 

already have been performed, as well as the Ricci scalar 

R = g af3 R Q/ g. (2.76) 

Then the five Einstein equations read: 

Goo = |(^'T+^^)-|r, R =0 
Goi = if), T il>,R--jjT, T = 

G 11 = l(ip, 2 T +i;, 2 R )-±r, R =o 



g 22 = -t,rr -r, TT -\ (-0,1 -ip,T ) = 

1 

R 



G33 = 2 ( i:, TT -i;, RR il>,n ) + G 22 = 0. (2.77) 
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We recommend the reader to distinguish the Roman R, the Ricci scalar, from the italic R, 
one of the Einstein-Rosen coordinates. 

Straight away one remarks that the component Goo is identical to Gn. Obviously only 
one of them is indispensable to obtain the complete set of the field equations, say Goo- 
Just this equation yields one of the differential equations determining if) and T, namely 
(I2.79p . Then of course we set in G 33 the component G 2 2 equal zero, so to get the Bessel 
differential equation (12.781) . From the component G i the last Einstein equation (I2.80p is 
read out. The result is a compact set of the Einstein field equations: 

^P,TT-^,RR-^,R= 

T, R =\R 

r, T = r if>, T ip, R 

This set of equations has a familiar structure: 

• Equation (12. 78j) looks exactly like the usual wave equation for the cylindrically sym- 
metric massless scalar field ip(T, R) propagating on a Minkowskian spacetime back- 
ground. 

• The succeeding equation ( 12.791) describes the energy density of the massless scalar 
field Tp in the cylindrical coordinates. 

• At last, through equation (I2.80p the radial energy current density is determined. 

As explained in the beginning of this chapter, the foliated (3+l)-dimensional general rela- 
tivity has been reduced to a (l+l)-dimensional model by the cylindrical symmetry. Further 
on, the model really is coupled to a massless scalar field, since such a field solves the Ein- 
stein equations written for the simplified cylindrical metric, namely the Einstein-Rosen 
metric. 

Last but not least we wish to notice that the set of the Einstein equations, reduced by 
the symmetry property of the Einstein-Rosen metric, can be deduced in a easier and faster 
manner. The whole calculation can be slashed to the task of complying 

R a(3 = 0. (2.81) 

The proof that the equations ( 12.711) are equivalent to the vanishing Ricci tensor is very 
short. The main step is to multiply (12.711) from the left with the Einstein- Rosen metric 
Oaf - One gets G = — R. This result is inserted in (I2.7ip wherefrom the following is given: 

— G a/ 3 = R a/ 3. (2.82) 

In our case of a vacuum spacetime the assertion R Qj g = follows trivially, as G a p = 0. 

Again the performed Ricci tensor for the Einstein- Rosen metric in Appendix C is taken 
for eliciting the set of the Einstein equations. From R 2 2 = as also from R 33 = we 



(2.78) 

(2.79) 
(2.80) 
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get the wave equation for the massless scalar field (12.781) . both equation are equivalent. 
Adding R 00 and Rn together and setting the sum equal to zero, (I2.79P is obtained. Finally, 
Roi = Rio = yields the energy current density (I2.80p of the scalar field. The main 
simplification in computing the Einstein equations with the second method is of course the 
fact that no Ricci scalar is needed. 



Chapter 3 

Einstein-Rosen Wave 



3.1 Solution of the Wave Equation 

In the present chapter we calculate the parameters if)(T,R) and T(T,R), which occur in 
the Einstein-Rosen metric and so determine the Einstein-Rosen wave. The cylindrically 
symmetric wave equation (12.78P is solved by a band of solutions. The parameter i(}(T, R) 
is a superposition of all allowed solutions. To keep track of, we note the fundamental wave 
equation again: 

The equation is obviously linear. The general solution of this differential equation is the 
cylindrically symmetric massless scalar field if}(T, R) propagating on a Minkowskian back- 
ground. The function ip(T, R) is the sum over all modes. A mode is an oscillation with a 
particular frequency u — k — p, (c = K — 1), k the wave number and p the momentum of 
the mode. As the frequency is a continuous parameter (k G Kg") the sum corresponds to an 
integration over all modes, i.e. over all k G K.q~, or equivalently over all momenta, p G WLq. 
Thus, the ansatz is the well-known Fourier decomposition by frequencies. The method also 
is known under the name "separation of variables" . Thus we arrange the following ansatz: 

1> h (T,R)=<Pk(F)XkW (3-2) 

This ansatz is put in the differential equation (13 .ip and thereupon we divide by the product 
tp k {T) Xk(R), in order to obtain 

fk,TT _ Xk, RR _|_ J_ Xk,R ^ 2^ 

X k R X k 

The ansatz (I3.2p presumes the functions (p k , \k t° depend only on T and R respectively. 
This fact constrains the right and the left side of the equation (I3.3P to be equal to a constant 
k 2 with k G Mq*. Otherwise the treated differential expression does not hold anymore. Then 
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we set 

-/,■- (3.4) 
-k 2 , k E Rq. (3.5) 



Xk, RR J_ Xk,R 

x k R x k 



The first relation is the ordinary free harmonic oscillator differential equation with oscilla- 
tion frequency k. This equation is solved by a linear combination of both solutions 

<4 1} Cn = A(fc)e* fcT 

V f(T) = B(k)e~ ikT . (3.6) 

So we find for the general solution of the time part of (13.1 j) the announced overlay: 

<p k (T) = A(k) e ikT + B(fc) e~ ikT . (3.7) 

The functions A(k) and B(fc) are in relationship with each other. As we know, gravitational 
waves, like electromagnetic waves, are real valued waves. From this it follows that the 
solution of the differential equation (13. ip has to be real. The function ifi(T, R) is real if and 
only if the following holds: 

B(fc) = A*(fc), (3.8) 

whereas A*(k) is the complex conjugate of A(k). Thereupon the solution to the time part 
of the wave equation reads: 

<p k (T) = A(fc) e tkT + A*(k) e~ lkT . (3.9) 

Subsequently we consecrate ourselves to the spatial part (13. 5p of the wave equation. 
Rewriting the equation into the form 

R 2 Xk,RR + RXk,R + k 2 R 2 X k = , (3.10) 
we notice the evident similarity to the Bessel differential equation 

d 2 d 

z 2 — w(z) + z — w(z) + (z 2 - n 2 ) w(z) = 0. (3.11) 
dz 2 dz 

In addition we set in (13.111) n = and in (I3.10p s = kR, s G M. The spatial differential 
equation is indeed the Bessel differential equation with index n — 0. 

s2 lR> Xk + S lR Xk + s2xk = °' (3 ' 12) 

with s = s(k; R). The Bessel function of first order solve the equations with n = 0. These 
are (see [5]): J (kR), Y (kR), H$f\kR), E { ~\kR), where tt^\kR) = J (kR) ± iY (kR). 
For the solution ip(T,R) the following boundary condition is given [3]: 

lim Ri), T ip, R = (3.13) 

H — ^0 
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Since the Weber function Y (kR) contains the logarithm function depending just on R, all 
expressions involving the considered function are invalid solutions, as the limes in (I3.13P 
does not exist. For this reason the cylindrically symmetric Bessel function Jo(fci?) is the 
last possibility to solve the spatial differential equation. For the further considerations 
we choose the representation for the Bessel function J (kR) given in appendix A, namely 
(IA.190 . but with z = kR and n = 0. Setting sin tp = — cosfi with ■& = ip + ~ and applying 
what is proved in (1A.12[) . we get for J (kR): 

j Q (kR) = — [ dye ikRcostp (3.14) 
2vr Jo 

Inserting J (kR) in (13.131) we obtain 

hm R cp k>T (T) J (LR) ip h (T) ^ MkR) = 0. (3.15) 

Thus, the Bessel function J (kR) solves the spatial differential equation and the resulting 
wave function if) k (T,R) the boundary condition too. 

ip k (T, R) = h{kR) [A{k) e ikT + A*(A;)e" ifcT ] . (3.16) 

The above function is the solution for a specific mode, namely the one with frequency 
uj = k. To obtain the general solution for the wave equation (13. 1ft . we have to sum together 
all modes, thus to integrate over k, accordingly to the Fourier theorem. Therewith the 
general solution of the wave equation for the cylindrically symmetric massless scalar field 
propagating on a Minkowskian spacetime is: 

POO 

tJ}(T, R)= dk Jo(kR) [A(Jfe) e ikT + A*(k) e~ ikT ] . (3.17) 
Jo 

If we compare the solution (I3.17P with the Hankel transformation ()A.22j) in appendix A 
for the case n = 0, we remark a close similitude. In fact the function ( 13. 17ft can be written 
as the Hankel transformation of ^y?fc(T), that is to say 



POO 

ip(T,R) = dkk J (kR) 
Jo 



A(k) JkT , A*(fc) „ 

-ikT 

k + k 



(3.18) 



Because of the close relation of ip(T, R) to the Hankel transformation, we define newly the 
expression -r(p k (T) . We set 

4>(k,T) = ]-<p k (T) (3.19) 



with the Hankel coefficients 



a(k) = a * {k) = ^1. (3.20) 
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Therewith the scalar field can be written as the Hankel transformation of i/j(k,T): 

POO 

ip(T,R) = / dk k J (kR) i){k,T) 
Jo 

POO 

= dkk J (kR) [a{k) e ikT + a*(k) e~ ikT ] . (3.21) 
Jo 

In order to be sure that the Hankel transformation of ip(T, R) is a well defined function, 
it is essential to give a look to the property the Hankel coefficients a(k) and a*(k) have 
to fulfill. As we can learn from appendix A, the two dimensional Fourier transformation 
with cylindrical symmetry is just the Hankel transformation - we can imagine that the 
Hankel transformation of ip(T,R) exists only if Tp(k, T), and in detail a(k) = ^jr— and 
a*(k) = A ^ , is integrable. In addition to this requirement, we would like that the 
function we get through the transformation, namely i/j(T,R), is again integrable over the 
radius R. This would assure the possibility to transform back i/)(T, R) to the functions 
a(k) and a*(k). This means mathematically expressed: 

7i : if)(R] T) — > ip(k; T) (Hankel transformation) (3.22) 

'HT 1 : "0(T; A;) — > tfj(R; T) (Inverse transformation) (3.23) 

and thus the identity 

H" 1 [H [ip(R; T)]] = H- 1 Up{T; k)] = ip(T, R). (3.24) 

In case of the Fourier transformation there is a function space S (Schwartz space), which 
guarantees the existence of the inverse map of a Fourier transformed function, if the treated 
function is element of S. This is ensured by the following theorem: 



Fourier inversion theorem: 

The Fourier transform is a bicontinuous bijection from 5(IR n ) onto 5(]R n ). Its inverse map 
is the inverse Fourier transform, i.e. (f)" = f = (f)~. The Fourier transform and the inverse 
map are denoted by f and f respectively [6]. 

Since the Hankel transformation is nothing but the two dimensional Fourier transformation 
with cylindrical symmetry it is assumed that an analogous theorem holds, which presup- 
poses the Hankel coefficients to be element of the Schwartz space in order to assure the 
identity (I3.24p . So, if it is set A(/c) and A*(k) to be element of the Schwartz space, then the 
functions a(k) and a*(k) are too, since the functions in S fall off faster than every power 
of their argument increasing to infinity. So we managed to assure the integrability of the 
wave function ip(T, R) and, of course, of its Hankel transformed function ijj(k,T) also. 



3.2 Finiteness of the Energy T(T, R) 

The integrability property of the wave function ip(T, R) is a vital condition for the function 
r(T, R). As it is shown in the appendix B, it turns out that the energy of the scalar field, 



3.2 Finiteness of the Energy F(T,R) 
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contained in a disk Az = 1, is given (up to a factor 2n) by this same function 

T(T,R)= I dR'T, R/ (T,R'). (3.25) 
Jo 

It is now important to prove that the energy remains finite for the radius R at infinit}0. 

POO 

T(T, oo) = I / dRR ) • (3-26) 

Jo 

Obviously, the properties of ip(T, R) play a main role for the finiteness of the energy. We 
first start with treating a single mode of ip(T, R) with the specific frequency k. 

MT, R) = Jo{kR) [A(fc) e ikT + A*{k)e~ ikT ] . (3.27) 

As we will integrate over R, we concentrate for the time being only on the f?-depending 
part 

MR) = MkR). (3.28) 
The plain wave is not normalizable as is generally knowro 

f°° 1 
WMRW = / dRRJ (kR)J ( y k'R) = -5( y k-k'). (3.29) 
Jo « 

In such cases a wave packet is introduced, which, for instant, could be formed by a Gauss 
profile A(k). This function is an element of the Schwartz space, and so falls off faster than 
every power of the argument k. 

POO 

ip(R)= dk A(k)3 (kR). (3.30) 
Jo 

Now the wave is normalizable, since the following integral is smaller than infinity due to 
the wave packet, which vanishes at infinity. 

POO 

\\tp(R)\\ 2 = dk \A(k)\ 2 < oo. (3.31) 
Jo 

Comparing f l 3 . 3 f) with ( 13.271) or directly with ( 13.171) we remark promptly that the wave 
function ( 13. 17ft is furnished already with its wave packet, and is therefore normalizable: 



||^(T,i?)|| 2 = / dRRi/; 2 (T,R) (3.32) 
Jo 

POO 

= dkk [a(k) e ikT + a*(k) e~ ikT ] 2 . (3.33) 
Jo 



1 The next equation is given by the expression (|2.79p for the energy density of the massless scalar field 

2 For eq. (|3.29j) : The orthonormality relation is presented in appendix A. Replace x with R in equation 
(TOST) . 
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So we can verify the finiteness of the energy T(T, R) if we integrate over the whole space 
at the fixed time T. For the derivatives of the wave function with respect to T and R we 
get 

POO 

ip, T = i dk k 2 J (kR) </?(_)(&; T) (3.34) 
Jo 

POO 

ip,R = - dk k 2 Ji(kR) <P(+)(k; T), (3.35) 
J o 



where we define 



V? (T) = a(k) e ikT =F a*(k) e~ ikT . (3.36) 
If we square both derivatives and integrate them over the radius R, we obtain: 

POO POO 

/ dR R ^ 2 , T (T, R) = - dk k 3 ip 2 ( ^(k;T) (3.37) 
Jo Jo 

poo POO 

/ dR Rtp 2 , R (T,R) = / dkk 3 if 2 {+) (k;T), (3.38) 
Jo Jo 

where we have used the orthonormality relation for the Bessel functional 

POO 1 

/ dR R J n {kR) J n (k'R) = —8{k- k ), Wn e Z. (3.39) 
Jo k 

Both integrals are put now in (13.261) . so to get the finite quantity 

POO 

T(T,oo) = 2 dkk 3 a(k) a*(k) 
Jo 

POO 

= 2 dkk 3 \a(k)\ 2 
Jo 

dkk\A{k)\ 2 , A(k),a(k) 6 S. (3.40) 

QED 

We conclude the present chapter affirming that the elicited wave function ip(T, R) is well- 
defined, since for the two cases, R — ► and R — > oo, the function remains bounded 
providing at the same time its normalizibility and a finite energy of the field. 



3 See (|A.26[) , for the case x = R. 



Chapter 4 
ADM Formalism 



In the previous chapters we have presented how the Einstein equations can be solved if the 
symmetries of the metric are known. In our case it is about the Einstein-Rosen spacetime 
and so the Einstein equations have been solved for the Einstein- Rosen metric (12.321) . The 
solution to the differential equations of second order is then, as we have seen, the massless 
scalar field In general the Einstein equations 

Ra/3~|^R=0 (4.1) 

are obtained by varying the following action with respect to the metric g a/3 , and discarding 
the resulting divergence^: 

S[g a p(x)\ = ^ j d i xC= ^ J ^Mx)R(x) d A x, (4.2) 

where £ stands for the Lagrangian density and g(x) for the determinant of the metric g a/3 
at the spacetime point x. In the Einstein equations above, R a p denotes the Ricci tensor, 
R the Ricci scalar and finally G the universal gravitational constant. In the early sixties 
of the last century R. Arnowitt, S. Deser and C.W. Misner, short ADM, developed the 
so called canonical formalism for general relativity [8]. A reason for constructing such a 
formalism to describe dynamic in general relativity is that an action containing canonical 
conjugate variables is a possible basis to develop a quantum theory of gravity by canonical 
methods. 

In order to introduce canonical variables, a frame is needed, where they can be related 
on. To this end, spacetime, in the ADM formalism, is split into three spatial dimensions 
and one time dimension, in literature shortly noted by (3+1). One can now attempt to 
imagine spacetime as a composite of hypersurfaces (3-space surfaces), which are fixed at 
arbitrary points in time. So a hypersurface Q tQ could be fixed at the time t = t , just as the 
hypersurfaces Q tk at an arbitrary time t = tk- Here we would like to consider two specific 
hypersurfaces Q tQ and Q tl to different times to < h- In the following we try to give an 
idea how points in the first hypersurface Qt are connected with the corresponding points 



^ee [7]. 
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in the later hypersurface Q tl . In the ADM formalism two functions are brought in, which 
measure the translation in spacetime. The lapse of the proper time between the 'lower' Q to 
and the 'upper' Q tl hypersurface is measured by the so called lapse function N(£, xi, x 2 , £3). 
On the other hand the shift function N*(i, xi, X2, £3) measures the spatial shift, i.e. this 
function specifies the position in the 'upper' hypersurface whereto the point x % of the 'lower' 
hypersurface have to be placed. The infinitesimal translation in spacetime is given by the 
interval ds 2 between the point x a = (t, x l ) and the point x 13 + dx 13 = (t + dt, x % + dx l ) in 
the 'upper' hypersurface, namely by the relation 

ds 2 = g^dx 1 + Kdt)(dx k + N k dt) - (Ndt) 2 = g a/3 dx a dx p . (4.3) 

Thus, the 4-metric has been constructed out of the 3-metric 3 gik, the lapse function and 
the shift function: 

' N,N S -N 2 N fe \ , N 

n, < 44 > 

If total divergences are discarded from the Lagrangian density of the action (14. 2p . then the 
functional S can be brought into the form 

dt J d 3 x Ng 1/2 (K ik K ik - K 2 + R). (4.5) 

Here, R is again the Ricci scalar and g 1 ^ 2 the square root of the determinant of the 3-metric 
go-. K 2 is the squared trace of the tensor K^: K = K* = g^K^. The new second order 
tensor Kjfc is the extrinsic curvature of the spacelike hypersurface. The extrinsic curvature 
is calculated out of the 3-metric, the lapse and the shift function, according to the formula 

K ifc = 7^( N i|fc + N fc|i - 9ik,o), (4.6) 

where the stroke means the covariant derivative in 3-space. As an illustration to the 
extrinsic curvature, we cite [2]: 

The extrinsic curvature measures the fractional shrinkage and deformation of a figure lying 
in the spacelike hypersurface S that takes place when each point is carried forward a unit 
interval of proper time "normal" to the hypersurface out into the enveloping spacetime. No 
enveloping spacetime? No extrinsic curvature! 

Varying the action (14.51) with respect to (7^,0 we obtain the momentum 7r jfc , which is 
canonical conjugate to the 3-metric gik- 

7T ik = = -gV 2 (K ik - Kg lk ), (4.7) 

ogik,o 

where the Lagrange density C is the integrand of the action (14.51) and 16/tG = 1 is set. 
Therewith it is scheduled what, in the ADM formalism, the generalized coordinates and 
the belonging canonical momentum are, namely the 3-metric gik and the momentum ir lk . 
Just the introduction of the momentum 
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makes it possible to transform the second order partial differential equations of the Einstein 
equations to an equivalent set of first order differential equations. This set will then depend 
on the dynamic variables gi k and n %k . The action functional is then converted into the 
Hamiltonian form, which is more convenient for the developing of a quantum theory of 
gravity by canonical methods, since the Hamiltonian enters the action. 

S = JdtJ d^x^g^o-NH-N^), (4.9) 

with S = S[<7ifc, 7T lh , N, Nj]. The lapse N and the shift function Nj in (14.91) play the role 
of Lagrange multipliers, whereby the function TC is the so called super-Hamiltonian and 
the function Ti 1 the supermomentum. These quantities are expressed by the canonical 
conjugate variables and 7r* fc , 

n = 9 - 1/2 (A,-i, 2 )-A, (4.10) 

W = -27r ik , k -g il (2g jl>k -g jk7l )^ k . (4.11) 

Analogously to K in ( 14. 51) . 7r denotes the trace of ir lk and so ir 2 is nothing but the square 
of the trace of ir lk : ir = 7r| = gi k n . The commas in the right hand side of the equation 
of the supermomentum represent the common partial derivative with respect to the space 
coordinates k,l = 1,2,3. Now, if the action (14. 9 j) is varied with respect to the Lagrange 
multipliers N and Nj, then the constraints of the system are obtained: 

H — 0, (4.12) 
W = Q. (4.13) 

Varying whereas with respect to the dynamical variables gi k and 7r* fc , the set of canonical 
equations is yield, which is well known from the Hamilton formalism and which, apart from 
the constraints, replace the Einstein equations in the form (14. ip . 

SR ^ m (A \A\ 

with the Hamiltonian 

B = J**<m + «#<) (4,5) 

We can see the relation (14.151) by recalling the definition of the Hamilton function: 

H = p i g i -L, (4.16) 

where p l are the canonical conjugate momenta to the generalized coordinates In our 
case we have then, 

L = J d'x {-K lk g lkfi -NH- N^ 4 ) , (4.17) 

P l qi = J d 3 xir ik g ikfi , (4.18) 
which reproduce equation (14. 15j) if inserted in (I4.16p . 
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4.1 Reduced ADM Action for the Cylindrically 
Symmetric Spacetime 

We dedicate this section to the ADM action for the cylindrically symmetric spacetime. Due 
to the symmetries of the model the 4-metric simplifies drastically, as it is shown in chapter 
one, see (12.2.11) . (12.2.21) . The first simplification is that the cylindrically symmetric metric 
depends only on the coordinates t and r, since it is invariant with respect to changes in 
the coordinates ip and z. Secondly, the metric acquires diagonal form by specific conformal 
transformations □, a fact which simplifies the task of solving the Einstein equations. These 
properties of course have influence on all quantities of the ADM formalism. We start with 
the lapse and the shift functions, which are related through equation (14.41) to the 4-metric, 
in our case the cylindrically symmetric one. The only component of the cylindrically 
symmetric metric, which in one hand is non vanishing and in the other hand is in relation 
with the lapse and the shift functions, is (I2.3ip 

g 00 = -n 2 (t,r), (4.19) 

where in the following we use the coordinates t and r in the conformal transformed metric 
instead of i and f. This component is only dependent on the coordinate t and r, such that 
also the Lagrangian multipliers N(i, r) and N 4 (t, r) of the ADM action in the Hamiltonian 
form only depend on these coordinates. Further on two components of the shift function 
drop out due to the symmetries, namely the ones which measures the changes in the <p and 
z coordinates. So, only the lapse function N(t, r) and one shift function N 1 (t,r), namely 
that one which measures translations in the radius direction, are left. Also the 3-metric 
depends only on the time and radius coordinates and as part of the conformal transformed 
metric, it is in diagonal form, see (I2.3ip . Since the canonical conjugate momenta are related 
to the 3-metric by the equation 

ogik,o 

one can imagine that the momenta inherit the properties of the 3-metric. In [3] the sym- 
metry reduction of the momenta is presented at length, so that we note only the results. 
The non vanishing components of the momentum tensor ir tk are the following, with 1, 2, 3 
labeling r, <p, z: 

7r n (t,r), 7r 22 (t,r), 7r 33 (t,r). (4.21) 

As expected, there remain only the components left, which are canonical conjugate to the 
corresponding non vanishing components of the 3-metric, 

gn(t,r), g 22 (t,r), g 33 (t,r). (4.22) 

If in the following discussion we take the metric (12.391) suggested in [3] as an example for 
a cylindrically symmetric spacetime, we see that the metric is given by the quantities 7, R 



'■see (pzm arm . 
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and The equations of motion, i.e the Einstein equations in the ADM form will finally 
be differential equations of those variables. So it is throughout correct to look at these 
quantities as generalized coordinates and to note the canonical conjugate momenta with 
respect to those generalized coordinates 

7 r 11 (*,r)=7r 7 e^-' y J vr 22 (t, r) = \RT 1 -k r 7r 33 (t, r) = (tt 7 + \Rn R + ir^e^, (4.23) 

where our calculation yields for the component TT 22 a factor R 1 instead of the factor R 
noted in the paper [3]. Thus, it has been achieved to bring the expression tt gik,o of the 
action (14. 9B in the following canonical form: 

n ik 9ik,o = tt 7 7 + n R R + n^ip, (4.24) 

where this form holds exclusively for the cylindrically symmetric spacetime and especially 
for the metric (12.391) depending on 7, R and ip. The dots in equation (I4.24I) denote the 
partial derivative d t . If (I4.23j) is solved for the momenta conjugate to the generalized 
coordinates 7, R and -0 then we get 

e^V 1 , (4.25) 
2 J Re" , V 2 , (4.26) 

^33 ^ _ n 22 R 2 e -i, _ ^11 e7 -^_ (4 27) 

Now it is straight forward to calculate the super-Hamiltonian and the supermomentum. 
To this end one inserts in ( I4.10p for 7r lfc the three non vanishing momenta TT 11 , 7T 22 and 7r 
and for the non vanishing components of the 3-metric, 

ffii = e^, g 22 = R 2 e~^ g 33 = e^. (4.28) 

The same is inserted in (14.111) for the computation of the supermomenta. To simplify 
the reduced ADM action anymore, the super-Hamilton and the lapse function as also the 
supermomentum and the shift function are rescaled respectively with the factors e 1//2< - 7- ^ 



and e^ 7 . The rescaled constraints read 

H = -7r 7 7r B + \R- l *l + 2R" - iR' + \R^'\ (4.29) 

Hi = -2^ + + ii'TTfl + VV 7 , (4.30) 

n 2 = n 3 = o, (4.3i) 



where the primes stand for the derivative with respect to the radius coordinate r. It is not 
surprising to see only one supermomentum, which does not vanish. We should take into 
account, there is only one non vanishing shift function, namely Ni(t, r). The variation of 
the ADM action in Hamiltonian form with respect to the rescaled Lagrange multipliers 
yields only the following two constraints, 



7T 7 = 
TtR = 
7Tw, = 



n = 0, 

fix = 0. 



(4.32) 
(4.33) 
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Inserting (14.241) and the rescaled constraints in the expression (j4.9p for the ADM action, we 
get finally the symmetric reduced ADM action for the case of the cylindrically symmetric 
spacetime described by the metric (12.391) . 

/oo roo 
dt dr (tt 7 7 + 7r fl fl + 7ty^ - Nft - N Hi). (4.34) 
■oo JO 



4.2 Reduced ADM Action for the Einstein- Rosen 
Wave 

After bringing the general cylindrically symmetric metric (12.31) into diagonal form by a 
conformal and a coordinate transformation into the inertial system, a further simplifica- 
tion of the diagonal metric was performed. By another conformal transformation of the 
coordinates^, the Einstein-Rosen coordinates were introduced. Exactly this step has to 
be repeated, i.e. to introduce the Einstein-Rosen coordinates in the canonical formalism 
such that the structure of the super-Hamilton and the supermomentum get more simpli- 
fied. The Einstein-Rosen radius R is already present in the canonical formalism, so that it 
remains only to enter the Einstein-Rosen time T. To this end the function T is defined by 
the extrinsic curvatur^] and is converted to a canonical coordinate by a canonical transfor- 
mation. The canonical coordinate T is then identified with the momentum 7r 7 , according 
to 

tt 7 = — T', (4.35) 

where the prime denotes the derivative with respect to r. The canonical transformation 
also allows to identify —7' in (14.291) and (14.301) with the canonically conjugate momentum 
of the Einstein-Rosen time, 

ttt = — V- (4.36) 

In section 12.3.11 we have seen that the introduction of the Einstein-Rosen coordinates 
changed 7 to the 'energy' function T of the massless scalar field following the relation 
(I2.43p . Through this transformation also the momentum ttt changes according to (14.361) . 
The result is the introduction of a new momentum Hp, which forces the implementation 
of a new momentum canonically conjugate to R, namely IT#. The relations between 
and I1 T and between tt r and 11^ are taken from [3]: 

IT T = 7r T +[ln(i?' 2 + T' 2 )]', (4.37) 



ln 



R' + r 

R'-T' 



-1 1 



(4.38) 



So the new set of the canonical variables for the reduced action is: (T, ITt), (R,Hr), 
(ip, 7ty). The canonical momentum is unaffected by the canonical transformation, since 



3 See section ()2.3[) and especially (|2.3.ip . 
4 See section VII in [3] 
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ip is untouched by it. Analogously ip did not change by the conformal transformation 

(fZEJ). 

Introducing the new canonical variables (T, Ht), (R,T1r) and (ip,7Tf) into the rescaled 
super- Hamilton (14.291) and the supermomentum (14.301) . we obtain: 

H = R'Ut + T'Ur + KR-^I + R^' 2 ), (4.39) 
Tit = T'Ut + R'Ur + iP'tt^. (4.40) 

Of course the constraints equations (14.321) and (14.331) remain unaffected by the reparametriza- 
tion, such that also for the reduced action in the Einstein-Rosen coordinates the constraints 
equation read 

H = 0, (4.41) 
Hi = 0. (4.42) 

The new canonical variables bring then the reduced ADM action in the subsequent sim- 
plified form, whereby S = S[T, U T , R, U R , if), n^}: 

dt / dr (n T T + Tl R R + -K^ip - N7Y - N Hi). (4.43) 

■oo JO 

By varying the action (I4.43P with respect to N and N 1 the constraint equations (14.411) and 
(I4.42p are derived. 



Chapter 5 

Observables of the Einstein-Rosen 
Wave 

In this chapter we dedicate us to the main topic of this work. It is presented how a set of 
phase space functions can be derived from the dynamical variables, namely the scalar field 
ijj and the canonically conjugate momentum tt^. The phase space functions turn out to be 
observables of the Einstein-Rosen waves. The condition the phase space functions have to 
fulfill in order to be considered as observables will be discussed in chapter six. 

Before we can look after the phase space functions, it is necessary to express the canon- 
ically conjugate momentum ir^ as a functional of the scalar field ip ( 13. 17ft . 



5.1 Canonical Momentum tt^ 

The canonically conjugate momentum is derived following the method presented in 
Torre's work [4]. In the canonical formalism the dynamic of a field is described by the 
Hamilton equation 

<A = {^H}, (5.1) 

where the dot denotes the time derivative dt- The scalar field if) is well known, being the 
solution to the wave equation (12.781) . The solution, we have found in chapter two, is 

POO 

tp(T, R)= dk 3 (kR) [A(k) e ikT + A*(fc) e~ lkT ] . (5.2) 
Jo 



The Hamilton function 



H= / dr(NH + N^i) (5.3) 



is given by the integral over the linear combination of the super-Hamiltonian 7i and the 
supermomentum 7i.i, which were obtained for the reduced action of the Einstein- Rosen 
spacetime. 

H = R'Ut + T'Ur + KR-'ttI + R^ 2 ), (5.4) 

n x = m T + #n K + v>V (5.5) 
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The coordinates T(r) and R(r) are part of an embedding X a : £ — > M 4 , S being the 
hypersurface at t — const., which maps the hypersurface S into the flat spacetime. The 
embeddings X a (r) = (T(r),R(r),tp,z) are arbitrary cylindrically symmetric slices on the 
interior of the spacetime and approach the hypersurface T = for r — > oo, where T is the 
Minkowskian time. The coordinates T(r) and R(r) do not depend on tp and z due to the 
cylindrically symmetric model. 

If the Hamilton equation (15.11) is solved then the canonically conjugate momentum 
is obtained as a functional of the scalar field ip. In the following prearrangements are 
presented, which allow to solve the Hamilton equation in a facilitated manner. It is then 
possible to write the time derivative, denoted by a dot in (15. ip . by a functional derivative 
of the embedding variables X a (r) = (T(r), R(r), tp, z). The definition of the unit vector on 
a given hypersurface is defined by 



x« a n a = o, 



aB 

9 p n a nf3 



(5.6) 
(5.7) 



where the index a runs over (1, 2, 3) labeling (R, <p, z). The metric g a p is the Minkowskian 
metric written for the cylindrically symmetric coordinates. The diagonal components are 
( — 1, 1, jR 2 , 1), while all other components vanish. The matrix X a a reads 

( T \ 
R' 
1 

\ o o i / 

From the condition (15.61) the two possible normal covectors on the hypersurface S are 
obtained: 



X 



,a 



(5.1 



n 



n. 



(i) 

a 
(2) 



(-R', T', 0, 0), 

(R',-r, 0,0). 



(5.9) 
(5.10) 



For the further calculation we choose without restricting generality the first one, which 
points out of the hypersurface in the negative time direction. The normal covector is 
normalized on —1 by the second condition (15.71) . 



n 

An a n 



A 



OR', r, 0,0), 

A(-R' 2 + T' 2 ) = 
1 



1, 



For convenience we assign the normalizing factor A to the covector, so to get in the further 
calculation a more compact form for the momentum 7ty(r), 

1 



n. 



rr 



R' 2 - T' 2 
(R',T, 0,0) 



-R',T', 0,0), 



(5.12) 
(5.13) 
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Hence, the induced metric on the hypersurface is given by the embedding variables, 

g ab = X« a9a pX\. (5.14) 

The components in the diagonal of the matrix g a b read 

(R' 2 -T ,2 ,R\l), (5.15) 

all other components are zero. The projection operator X a a , which takes a spacetime 
vector into a vector in the hypersurface, is expressed by the induced metric, 



X\ 



a&X , 



9a(3 



b g 



I 


— T" 





o \ 


R' 2 -T' 2 




R! 










R'2_ T ,2 







1 





\ 








1 / 



(5.16) 



If the constraint functions (15.41) and (15.51) are taken together to a covector then it is possible 
to express the Hamilton functions as follows: 



H 



drN H a , 



(5.17) 



where a = 0, 1,2,3. Of course 7Y 2 = = 7i 3 due to the symmetries, as we have learned 
in section (14.11) . In this notation the constraint equations read: 7i a = 0. Then also the 
Einstein-Rosen momenta and in H5.4I) and H 5 . 5 f) respectively canonically conjugate 
to the embeddings T(r) and R(r) are taken together to a covector, 



n Q = (n T ,n R ,o,o). 

Thus the Hamilton function in (15. 17ft is newly expressed by, 

H a = U a -g- 1 / 2 n a h + X a a h a , 

h = \ (vrj + R 2 tp' 2 ) , 
h a = (tt^',0,0). 



(5.18) 



(5.19) 

(5.20) 
(5.21) 



The constraint function in this form allows to solve the Hamilton equation in a easier way. 

(5.22) 



V>(*,r) = |^(*,r),jf dr'N (r')H a (r') j . 

In order to compute the Poisson bracket, it is necessary to choose a lapse and a shift 
function. In our case the fixing of the N is not needed, as the introduction of embeddings 
selects a precise folium (an instant of time) i.e. t = const., what just corresponds to the 
choice of N . In fact, only the r-dependence of the quantities will be noted, for example 
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ip(r) instead of ip(T(t,r), R(t,r)). The choice of the embeddings permits to replace the 
time derivative by the functional derivative. 



oo / j \ r poo 



d r°° . / d 



, dr —X a (r) t--—- = dr N (r) (5.23) 
dt J \dt KJ J5X«(r) J v ! 5X«{rY v ; 

where it has been utilized: 4X a (t, r) = N . If the operator fl5.23f) is applied on the time 
derivative of the wave function, one gets 

far) = J™ dr' N>') = dr' N°(r') {^(r), Hjf)} . (5.24) 

The right hand side is obtained by simplifying the same side of the Hamilton equation 
(I5.22p . considering that Poisson brackets are linear in both arguments, so 

{POO ~\ POO 

ip(r),J dr'^ a {r')H Q {r')\ = J dr' N°(r'){^(r), H a (r')}. (5.25) 

Comparing the integrands in (15.241) . the following has to hold obviously. 

^)={#)-^M}. (5-26) 

If in the above equation the Poisson bracket is solved then the conditional equation for the 
canonical momentum 7r^,(r) comes out. 



Mr),W a (r')}= / dr 



oo 

II 



5ip{r) 8H a {r') 54){r) 5H a {r') 



8ij}{r") 5n f (r") 5^{r") 5^j(r") 
dr" 5(r — r") 8{r — r') [—g~ 1 ^ 2 n a 7T^{r') + X^tp'ij- 1 )] 
= -g' 1/2 n a TT^(r') 5(r - r') + X\ ip'(r') S(r - r'). (5.27) 
The canonical equation ip = {ip, H} has been reduced to 

<)( ' (r) S ( r _ /) [-g-^n^ir') + X>'(r')] (5.28) 



5X a (r' 



To get the momentum, we set this form of the canonical equation in ( 15. 24ft . remembering 
g^ 1 / 2 = R^ 1 , g is the determinant of the Minkowskian metric in cylindrical coordinates. 

n a {r) = n « (r ) X i^ {r) _ R -l n a na %4r y (5 2g) 
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Since the smear function N (r) can be arbitrarily chosen, it is allowed to substitute it with 
the unit vector n a . As n a n a = — 1 and n a = 0, the formula for the momentum finally 
follows from (15T29D : 

^(r) = R(r) [R r Tp?(T(r),R(r))+T r ^ R (T(r), R(r))]. (5.30) 

The canonical momentum -rr^ is then given, as suited, by a functional of the massless scalar 
field tp. Inserting the derivatives with respect to the embeddings the momentum is then 
written as a functional of the mode expansion A and A* contained in the field ip. 



5.2 Cauchy Data 

The cylindrically symmetric wave equation ( 12.781) as also equivalently the Hamilton equa- 
tion ( 15. ip define the initial value problem for the scalar field ip. If initial data are given, 
also known as Cauchy data, the dynamics of the field ip(T, R) is uniquely described by the 
differential equation (12.781) . If the problem is solved in the canonical formalism then an 
initial value of the field ip(T, R) and of the momentum ir^(T, R) must be given to uniquely 
fix the dynamics of the field by the conditional equation of the field and momentum. The 
field ip(T, R) is known as solution to the wave equation ( 12.781) . see for instance (15. 2D . The 
momentum tt^(T,R) is obtained by inserting ( 13.341) and (I3.35f) in equation (15.301) . Then 
the momentum 7i^(T, R) reads: 

POO 

Tc^p (T, R)=iR dkk J (kR) [A(k)e ikT - A*(k)e~ ikT ] R, r 
Jo 

POO 

-R dkk h(kR) [A(k)e ikT + A*(k)e~ ikT ] T, r . (5.31) 
Jo 

Now we affiliate the Cauchy data ipo(r) an d n^ (r), which are calculated for T{r) = and 
R(r) = r: 

POO 

Mr) = ^{T, R)\t=o,r=t = / dkJ (kr)[A(k)+A*(k)}, (5.32) 

Jo 

POO 

^„(r) = ir^T, R)\ T=0!R=r = ir / dk k J (fcr) [A(fc) - A*{k)} . (5.33) 



In order to enter the initial data ipo{r) and 7ty (r) in the general solution ip{T,R) and 
7r^,(T, R), ip^ir) is adapted to 7r^ (r). Subsequently the sum of both is Hankel transformed. 



ikrip Q (r) = ir dk k J (kr) [A(k) + A*(k)} 
Jo 

POO 

7r^, (r) + ikrip (r) = 2ir / dk k 3 (kr)A(k) 



o 



2iA(k) = / dr J (fcr) [^ (r) + ikrip Q (r)\ ■ ( 5 - 34 ) 
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Hence, the expansion coefficients A(k) and A* (A;) are given as functionals of the Cauchy 
data ipo( r ) an d 7ty ( r ) f° r T{r) = and R(r) = r. 



A(k) 
A*(k) 



dr J (Ar) [krip (r) - i7i> (r)] , 
dr J (Ax) [kri(j (r) +Z7r^ (r)] . 



(5.35) 
(5.36) 



By inserting the mode expansions (I5.35P and (15.361) in the general solution (15.21) and in the 
momentum ( 15.311) . we managed to express the canonical quantities tfj(T,R) and tv^(T,R) 
as quantities of the Cauchy data V'o( r ) and 7ty (r). 



ip [A(k),A*(k); T(r), R(r)) 



dkJ (kR) 



dr J (Ar) (krvp -iir^e 



ikT 



+ 



dr J (Ar) (kri/j +1^] 



-ikT 



(5.37) 



poo 

tt^ [A(k),A*(k); T(r), R(r)] = iR dkk J (kR) 

Jo 



dr J (At) (kr ^ - i 7T^ ) e 



ikT 



dr J (Ar) (kr i{j + iir^ ) e 



-ikT 



RlT 



-R dkk J^kR) 



dr J (Ar) (krip Q -itt^) e 



ikT 



+ | / dr J (kr) (krip + iTV^) e 



-ikT 



T 

1 V 



(5.38) 



Hence, the dynamics is uniquely given by the initial data 4>o(r) and vr^, (r) and so the 
initial value problem has been solved completely. 



5.3 Observables A(k) and A*(fc) 

In the final part of the present chapter we would like to deduce the phase space functions 
A(k) and A*(k), which turn out to be observables of the cylindrically symmetric gravi- 
tational waves. By deriving the field ip and its conjugate momentum as functional of 
the Cauchy data (I5.32f) and (15.331) . we also calculated the phase space functions A(k) and 
A* (A) for a specific hypersurface, namely for T(r) = and R(r) = r. Of course, our goal 
is a generalization of these mode expansions, that is to say for the phase space functions 
of an arbitrary hypersurface. To this end, we perform some Hankel transformations of the 
canonical variables x/j and n^, whose sum finally will lead to the desired set of observables. 



5.3 Observables A(fc) and A* (A:) 
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In the beginning we start to transform the momentum while multiplying it with a phase, 

POO 

i / dre ikT 3 (kR)n4r) = - e lkT [A(k)e lkT - A*{k)e~ lkT ] 
Jo 

! POO 

drR dkkJoik'RjJ^kR) 



x 



[A(k)e ikT + A* (k)e~ ikT ] T, r e ikT . (5.39) 



The field is treated in a similar way, in order to obtain the integral term in the right hand 
side of the equation above. 

poo poo poo 

i dre ikT "fc Ji(fcJ2)T, r ^(r) =i / drR dk k J (k' R)J 1 (kR) 
Jo Jo Jo 

x [A(k)e ikT + A*(k)e~ ikT ] T, r e ikT . (5.40) 

Further on ip(r) is transformed once more in order to get something related to the first 
part of the transformation (15.391) . 

POO 

/ dre ikT kRR,rJ (kR)ij(r) = e lkT [A(k)e lkT + A*(k)e~ ikT ] . (5.41) 
Jo 

If the sum over the three transformations is taken, we easily obtain the observables. 

POO 

2A*(k)= dre ikT {ij(r)kR[iT, r J 1 (kR) + R, r J {kR)} + iJo{kR)7i^}. (5.42) 
Jo 

Thus, the observables on an arbitrary hypersurface are: 



A(k) 



POO 

\ / dre- ikT {ip{r)kR[R, r J {kR)-iT, r J 1 (kR))-iJ (kR)7r f },(5A3) 
Jo 

POO 

A*(k) = \ \ dre ikT {ijj(r)kR[R, r 3 (kR)+iT, r 3 1 (kR)} + i3 (kR)7i i ,}. (5.44) 
Jo 

These phase space functions yield just the mode expansions (15.351) and (15.361) if T(r) = 
and R(r) = r are set in the observables above. So (15.431) and (15.441) are really generaliza- 
tions of the phase space functions ( 15.351) and ( 15.361) . 



Chapter 6 

Algebra of the Observables 



For the developing of a quantum theory of the cylindrically symmetric gravitational waves, 
observables of the system are needed, which enter the theory as operators. So it is of 
crucial importance to investigate the gauge invariant property of the phase space functions 
proposed to be observables of the model. Then, if the derived phase space functions A 
and A* on an arbitrary hypersurface are suggested to be observables of the treated model, 
gauge invariance of these functions has to be assured. In the first part of this chapter we 
present the condition a phase space function has to fulfill in order to be considered as an 
observable. Then we will show that the phase space functions A and A* are observables. 
In the second part we will calculate the Poisson brackets of the observables in order to 
derive their algebra. 

6.1 Gauge Invariant Phase Space Function 

A classical observable is a phase space function, which is gauge invariant on the constraint 
surface [9]. First class constraints 7 a of a system are quantities, which are the generators 
of infinitesimal gauge transformations. A phase space function, which remains constant on 
the orbit in the constraint surface is gauge invariant. This is equivalent to the vanishing 
Poisson brackets of the phase space function F and the first class constraints -y a , 

{F, 7 a}~0. (6.1) 

The bracket is calculated on the whole phase space and afterwards its value is taken on 
the constraint surface. The weak equality symbol is to be understood in this way. In 
the following section an example for this kind of calculation is given. If (16.1 j) is fulfilled 
then we are allowed to consider the phase space function F as an observable of the system. 

6.1.1 Gauge Invariance of A(k) and A*(k) 

In order to find out whether the mode coefficients A(k) and A*(k) of the cylindrically 
symmetric scalar wave are gauge invariant, we have to assure that the Poisson bracket 
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of these functionals with the first class constraints vanishes on the constraint surface, 
according to the formula (16.11) . The first class constraints are the super-Hamiltonian (15.4)) 
and the supermomentum ( 15.5)) . Thus, A(k), and therewith also A* (A;), are observables if 
the following formula holds: 

{A(k),H a (r)}^0. (6.2) 

The constraint submanifold T is the surface where TC a = 0. Then, saying the Poisson 
bracket (16. 2p has to vanish on the constraint surface means that one evaluates firstly the 
bracket and secondly, if still necessary, Ti a = is set, so to investigate whether the bracket 
really vanishes. The relation (16.21) can also be rewritten equivalently as, 

{A(fc),W a (r)}| r ^ a=0 = 0. (6.3) 

As shown in [4] the Poisson bracket ( 16.21) vanish, such that the phase space function are 
observables. The calculation is performed for a scalar field propagating in curved spacetime. 
The cylindrically symmetric scalar field on a Minkowskian background is then a special 
case of the presented computation. 



6.2 Poisson Brackets 



In this section we are going to calculate the Poisson brackets of the phase space functions 
A(k) and A*(k). The motivation therefore is to elaborate a Poisson algebra of these 
observables. To this end we introduce the functions Q(r) and P(r) written in [4], where 
just for consistency we retain our normalization. 



Q(r) 
P(r) 



dk J (fcr) [A(fc) + A*(Jfe)] 



ir / dkkJo(kr) [A(k) - A*(k)] 



(6.4) 
(6.5) 



where Q(r) = ipa(r) 05. 32)) and P(r) = 7Ty, (r) ( 15. 33f) . As we know from section l5\2l the 
functions Q(r) and P(r) are related to A(k) and A*(k) by a Hankel transformation, which 
will be undone in the next step so to describe the phase space functions as functionals of 
Q(r) and P(r). 



A(k) 
A*(k) 



dr Jo(kr) [kr Q(r) — i P(r)] 
dr Jo{kr) [k r Q(r) + i P(r)] 



(6.6) 
(6.7) 



6.2.1 Algebra of A(k) and A* (AY) 

It is straight forward to compute the Poisson brackets for the observables A(k) and A*(k). 

~8A(k)5A*(k') 5A(k)5A*(k')~ 



{A(fc),A*(fc')} = / dr 
'o 



5Q(r) 5P(r) SP(r) 5Q(r) 



(6.i 



6.2 Poisson Brackets 
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5A(k) 



\kr J (£x), 



5A(k) 
5P(r) 



(6.9) 



<JA*(fc') 



Jo(fcV) 



— \kr J (k r). 



(6.10) 



<SP(r) 2 " uv " / ' 5Q(r) 

Inserting the functional derivatives in (16.81) . we obtain the coefficients A(fc) and A*(k) to 
be observables, which are canonical conjugate to each other. 

POO 

{A(Jfe), A*(k')} = dr[\ kr J (fcr)J (fcV) + \ k'r J (fcV) J„(Ax)] , 
Jo 

= lk\5{k-k') + {k'^5{k-k'), 



k 

S(k-k'). 



(6.11) 



{A(fc),A(fc')} = / dr [-{kr J (kr)J (k'r) + { k'r J (fcV) J (Ar)] 
Jo 

= -{k\8{k-k')+{k'y8{k-k'), 



0. 



(6.12) 



{A*(A;),A*(A;')} = / dr [{k'r J (kr)J (k'r) - \kr J (kr)J (k'r)] 
Jo 

= ik'y5{k-k')-{k\8{k-k'), 



k 



0. 



(6.13) 

Hence, the well-defined algebra for the phase space functions A(k) and A*(k) of the cylin- 
drical gravitational wave has been elicited: 

{A(fc), A*(fc')} = 1 8(k - k'), {A(k),A(k')} = 0, {A*(fc), A*{k')} = 0. (6.14) 
6.2.2 Algebra of Q(r) and P(r) 

Analogously, the Poisson brackets of Q(r) and P(r) can be calculated, so to obtain the 
Poisson algebra of these observables. 

'tfQ(r) (JP(r') <JQ(r) 5P(r') 



{Q(r),P(rO}= / dk 



SQ(r) 
5A(k) 



h{kr) 



_8A{k)5A*(k) 5A*(k)5A(k) 
8Q(r) 



SA*(k) 



h{kr) 



(6.15) 
(6.16) 



SP(r') 
5A{k) 



ikr' Jo(fcr'), 



5P(r') 
5A*(r) 



-i k r' Jo(kr'). 



(6.17) 
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From this it follows immediately the Poisson algebra for Q(r) and P(r): 

{Q(r),P(r J )} = -2i5{r-r r ), {Q(r), Q(r')} = 0, {P(r), P(r')} = 0. (6.18) 

It is not surprising to get the same Poisson algebra for A(fc) and A*(k) as well as for 
Q(r) and P(r) until normalization and sign. The functions Q(r) and P(r) are, from a 
mathematical point of view, just expressions for the Hankel transformed real and imaginary 
parts of the mode coefficients A(k) and A*(k). 



6.3 Quantum Theoretical Property 

The Poisson algebra ( 16. 14ft for the observables A(k) and A*(k) has an interesting feature 
if investigated from a quantum theoretical point of view. As known the Poisson brackets 
multiplied by the complex % correspond to the commutators of quantum theory. Multiplying 
the Poisson brackets in (16.141) with the complex i, the commutation relations for the creator 
and the annihilator operators of the scalar field - see for instant [10] - are obtained up to a 
normalizing factor a/2. Hence, the observables A(k) and A*(k) can just be identified with 
the creator and annihilator respectively. 

A -f A, 

A* -> A 1 ". (6.19) 

Viewing now the cylindrically symmetric scalar field as a wave function in quantum field 
theory and applying the operators ( 16. 19ft and ( 16.191) on the wave function, the Fock space 
of the system can be constructed, which spans the Hilbert space of the quantum states of 
the cylindrically symmetric scalar field. 



Chapter 7 
Conclusion 



By a conformal deformation and a coordinate transformation into the inertial system the 
general cylindrically symmetric metric is brought into diagonal form. The introduction of 
the Einstein-Rosen coordinates allows a further simplification of the metric. The symme- 
tries and the special form of the metric reduce the Einstein equations to a set of three 
differential equations. The master equation is solved by the cylindrically symmetric mass- 
less scalar field propagating on a Minkowskian background. 

The introduction of the reduced ADM action for the Einstein-Rosen wave, and there- 
with of the super-Hamilton and supermomentum gives the opportunity to derive the 
momentum, which is canonically conjugate to the scalar field, by solving the Hamilton 
equation of it. Combining the scalar field with its canonically conjugate momentum and 
performing several Hankel transformations it is possible to isolate the mode expansions 
contained in the function of the scalar field. Just these phase space functions turn out to 
be gauge invariant and hence observables of the model. This is shown by investigating the 
Poisson bracket of the phase space function and the first class constraints. As it is required 
for observables the Poisson bracket vanishes. 

In order to guarantee the phase space functions to be canonical conjugate to each 
other the Poisson brackets of them are computed. The result assure the expected property 
and allow the correspondence of the observables to the annihilator and creator operators 
in the formulation of the quantum theory for the cylindrically symmetric scalar field on a 
Minkowskian background. Acting on the scalar field, now considered as a wave function, 
the Fock space can be developed, spanning therewith the Hilbert space of the quantum 
states of the cylindrically symmetric scalar field. 



Appendix A 



Orthonormality Relation for the 
Bessel Function of the First Kind 



A.l One Dimensional Fourier Transformation 
in Cartesian Coordinates 

In the beginning we shall study the following integral: 

POO 

dk e -ikx e ikx> = dkel (x-x>)k_ (A.l) 



Unfortunately this integral is divergent! In order to obviate this trouble, a suitable regu- 
larisation is carried out. To this end a new measure is defined, which contains a proper 
damping, whereas e > 0. 

[dk] = (0(-Jfe) e tk + 6{k) e~ ek ) dk. 
Therewith the integral exists and can be calculated. 

/0 poo 
dk e ek e i(x-x')k + / dk e -ek e i(x-x')k 
-00 J 

1 1 

with u = x — x 



e — i(x — x') e + i(x — x') 
2e 



e 2 + u 2 



The next step is to perform the limes e — > to achieve the value of the initial integral 



/OO POO 
dke i(x-x')k = Hm / [dk] e^-^ 
■00 J-00 



—00 
2e 



lim = 2vr 5(x - x') (A.2) 



e^o e z + it 



J In [11] we find: lim e ^ ^^4^ = <5(x), for e > 0. 



48 



Chapter A. Orthonormality Relation for the Bessel Function of the First Kind 



If the distribution flA.2[) is applied on a test functional e S, so obviously the result is not 
normalized. To adjust this non-conformance it is sufficient to assess the integral flA.ip with 



j-. Now we are allowed to undertake the identity map of the Fourier transformation in 
one dimension. Let's define the transformation and its inverse respectively as follows: 

i(x) = — / dke ikx i{k) 



2tt 



oc 



/oo 
dx' e~ ikx ' i(x') 
-oo 

The identity map reads 



f( x ) = _L/ dx'e lkx l dke- ikx 'f(x'), (A. 31 



oo J — oo 



whereby the integration order has been changed, since the integral over k is regularizable, 
as shown above. So we obtain 



-I /*oo ^oo 

f(x) = — dx' dk e i( - x - x ' )k f(x') 

J —oo J —oo 



oo 



— / dx' 2u 5{x - x') i(x') = f(x). (A.4) 



2tt 

Finally, the identity map ( 1A.3I) is normalized on one by the factor Nd=i = 

A. 2 Two Dimensional Fourier Transformation 
in Cartesian Coordinates 

In order to understand better the switch from the transformation in Cartesian coordinates 
to the one in cylindrical coordinates, first we shall write down the relations in the first 
mentioned coordinates. 

We take two vectors from the two dimensional vector space Mr, namely X — (^1 5 ^2 

) and 

k = (ki, k 2 ). The Fourier transformation in d = 2 and its inverse are defined respectively 
as follows: 

f (^) = 77^ / d 2 ke* kx i(k) (A.5) 



oo 



f(fc) = / d 2 x'e- ikx " f(f). (A.6) 



2 iS is the Schwartz space. 



A. 3 Hankel Transformation 
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Applying the inverse on the transformation we get as is well known the identity mapE): 

-| /*oo -I roc 

f(f) = _ / d 2 ke 1 ^— / d 2 x' e" ife ' f(f) (A.7) 
2tt 2tt 

OO POO ^ 

d 2 x' / d 2 k e ^-^ k f(x') 



(27T) 

/OO /*00 1 /*Cxj 

-OO »/ — OO ^TT »/ — OO 



— OO 



(2tt)' 

f(f). 



1 r°° 

x — / dk 2 e l{x2 - x 'i )k2 \x[,x' 2 
dV (2tt) 2 5 2 (f-f)f(f) 



(A.8) 



In the separated representation of ( 1A.8I) we observe twice a one dimensional integra- 
tion over ki resp. k 2 of the form ( 1A.1I) . Belonging to ( 1A.2I) each of them are equal to 
27r 5(xi— x'j), i = 1,2. Consequently, we understand the reason why in ( 1A.5I) the normalizing 
factor in two dimensions A^ =2 = appears, analogously to the Fourier-transformation 
in one dimension. 



A. 3 Hankel Transformation 

We start now with noting the two components of both vectors x = (xi, X2) and k = (ki, k?) 
in cylindrical coordinates. That is 

£1=0; cos a, X2 = xsma ki = kcos/3, k 2 = k sin /3. 

For further need, we prepare a third vector from R 2 in the same coordinates: x! = (x[, x' 2 ) 
with x[ = x'cosa', x 2 = x'sina'. Thus, we get for the following products: 

kx = kxcos(a — (3), kx' = kx' cos(a' — (3). 

Please note that the scalar s k, x and x 1 G R are the absolute values of the homonymic 
vectors. Furthermore the determinants of the Jacobi matrix for the integral over x' and k 
give the measures dx' times x' and dk times k. A function g(x) = g(x%, x 2 ) and its Fourier 
transformed function g(x') = g{x' x , x 2 ) can be noted in polar coordinates as follows: 

g(x) = h(x, a) = f(x)e ma (A.9) 
g(x') = h{x', a') = f(x')e w . (A.10) 



3 Belonging to (j A.2|) the integral (|A. 1|) is regularizable. Thus, the integration order is arbitrary. 
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If we write down the Fourier transformation identity map in cylindrical coordinates, then 
it reads: 

f( x ) e ina = —— dkk d(3e lkxcos ^- a) dx'x' da' e~ ikx ' ^ ~ 0) e w i(x') 
(2tt) 2 J J_„ J J_ n 

(All) 

We treat now the last equation in order to obtain finally the identity flA.7|) containing the 



Bessel function J n (kx). For this purpose we prove that an integral of the form 

«2tt 

dip e ifc:ccos (^-W)) 



o 

is equal to 

r2n 



This is possible due to the symmetry property of the trigonometric functions if a whole 
period is taken. That is gladly our case. 

2tt f^O p2n 

dip e ikxcos (f-'Po) — / (fan gifescosC^-Vo) _j_ / fa e ikxcos(ip-ip ) 

(I JO Jtpo 

2it p2iv-\-<pQ 

d(f e ikxcos (<fi-<fio) _|_ / ^ e ikxcos(<p-ipo) 

ipo J 2n 

2-K+ipo 

dip e ikxcos (<fi- l P") 

2n 

d$tf kxcos *, (A. 12) 

o 

whereas •& = ip — ipo- 

QED 

Thereafter we return to the crucial equation ( 1A.11I) . We perform a first substitution of 
variables, namely 

/?-a=|+/3'. (A.13) 



The transformation map flA.lip changes to 



oo poo i"K—a 



{{x)e ina = / dx'x' / dkk dp'e~ ikxs ' m P' 

(2tt) j Jo Jo J-^-a 

x / da'e- ikx ' cos{a '~ a - p -^e ina '\x') (A.14) 



In doing so we have employed the properties of the cosine function 

cos((/?) = cos(— ip) (A. 15) 

7T 

cos(- + ip) = -sin(y?). (A. 16) 



A. 3 Hankel Transformation 
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In virtue of the periodicity of the trigonometric functions, the interval [— |tt — a, | — a] 
can be traced back to the original interval [— 7T, n], cp. (IA.12p . Then a second substitution 
is undertaken, namely 

a ' - a - = ■&, (A.17) 
whereby one gets for the transformation newl)Q: 

/*00 /*00 /»7T 

f(x)e inQ = e ina / da;' a;' / dkkj^ rf/ g' e - fcsin/3 ' +in/3 ' 



(2t) , 

^0 ./— 7T 



x ^ / dde- ikx ' s[nM ^~f(x'). (A.18) 



(™7T 

Also for this step the periodicity of the trigonometric functions has been of use, since 
one integrates at each case over a complete period. Further on the identity cos(i9 — | ) = 
cos(| — i?) = sin $ has been utilized, too. The integral representation of the Bessel function 
J n (z) is given in [12] by 

J n (z) = J_ / dtp S"*-*^ . (A.19) 

27T ./_,r 



By setting i? = ^ + 7r and making use of cos(<^ + 7r) = — cos and sin(y9 + 7t) = — sin 
it can be shown that the Bessel function 3 n (z) is a real function, i.e. J n (z) = ^n( z )- 
Comparing both integrals over (3' and d in ( 1A.18I) with the integral representation of the 
Bessel function J n (z), we remark that the integrals in (1A.18I) are just the Bessel functions. 
After latter cognitions the sums over (3' and d can be substituted with the Bessel function 
J n (kx) and J n (kx'), replacing z with kx respectively with kx' . Thus, we get: 

f(x) = / dx' x' dkk 3 n (kx)3 n (kx') f(x') = f(x) (A.20) 
Jo Jo 

If we are interested in the Fourier transformation on the dependence of the vector's x 
absolute value only, then the function f in (1A.18I) is no more subject to the angle a. 
Recapitulating we put down: 



g(\x\) = f(x) = / dx' x' dkk J n (kx)J n (kx ; ) f(x') = f(ar). (A.21) 



Analogously to the identity map in Cartesian coordinate ( 1A.7I) . the Fourier transformation 
and its inverse, described in the new variables, can be obtained comparing just the form 
in (1A.7|) and the gained form in cylindrical coordinates in ()A.20j) . The so-called Hankel 



oo 



transformations read: 

f(x) = I dkk J n (kx) £(k) 
Jo 

POO 

f(fc) = / dx' x' J n {kx') f(x'). (A.22) 







4 Thc integral succession can be switched arbitrary, as the regularisation for the one dimensional case 
in section A.l guarantees the existence of the integrals. 
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These equations are called Hankel transformation. They are the analogous expressions to 
the two dimensional Fourier transformation with cylindrical symmetry. The normalization 
factor for this transformation is also given in the last two equations, namely N cy i coor( i. = 1. 
So we deduce that the Bessel function defined in [12] already are normalized. Observing 
properly fjA.llj) . we remark that the initial factor -^p- has been absorbed in both Bessel 
functions J n (kx), J n (kx r ) belonging to the definition in [12]. 

As a main result of this appendix we want to write down explicitly the normaliza- 
tion factor for the orthonormality relation of the Bessel function of the first kind J n (kx). 
We can take the relation directly from the expression ( 1A.20I) : 

poo poo 

f(x) = / dx' x' dkk J n (kx)J n (kx') i(x') = f(x) (A.23) 



The orthonormality relation is just the integral over k, which obviously has to be equal to 
the Dirac 5-function multiplied with the normalizing factor. 

1 

dk k 3 n (kx)3 n (kx') = — S(x — x), (A. 24) 

x 

with the normalization factor 

Nx' = - r (A.25) 
For the Hankel transform f(k) there is of course the analogous orthonormality relation: 

dx x J n (kx)J n (k'x) = — S(k — k'), (A.26) 

k 

with the corresponding normalization factor 



(A.27) 



Appendix B 

The Energy F(T, R) 



The metric (12.321) is determined among other things also by the function T(T,R). This 
function is given indirectly by the gradient formed by the derivatives T,t (T, R) (12.801) and 
T, R (T, R) (12.791) . which describe the energy density of the massless scalar field ip(T, R) on 
a Minkowsikian spacetime background and respectively the radial energy current density. 

VT(T,R)=( 1> T ) (B.l) 



r. 



In order to get the function T(T, R) we have to integrate the gradient field (IB. II) . Starting 
from the point T(0,0) we would like to reach the arbitrary point T(T,R) by choosing an 
arbitrary path. However, at this place we make use of a fundamental property of gradient 
fields. This sort of fields are conservative. This property implicates that all path from 
the starting point T(0,0) to the arbitrary point T(T,R) are equivalent. So we are free to 
choose a particular path without restricting generality. 

r{T,R) = r(o,o) + f dr r> (t',o) + fan! r, R , (t,r'). (b.2) 

Jo Jo 

We prove now that 

T(T,R) = f dR' F,ri (T, R'), (B.3) 
Jo 

if we deal in the following with a pure radiative field without sources on the axis of sym- 
metry (z-axis). 

Proof: If we have no sources on the axis of symmetry, the spatial geometry must be lo- 
cally Euclidean on the mentioned axis. So the circumference of a small circle i?=const., 
T=const., z=const. is the 2tv multiple of its proper radius. For i?=const., T=const., 
2=const. the Einstein- Rosen line element 

g%K(T, R) x a = e r "^ {-dT 2 + dR 2 ) + R 2 dip 2 + e^ dz 2 (B.4) 
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reduces to the formula 

ds 2 = R 2 d<p 2 (B.5) 

^ ds 2 = Re~^ dtp. (B.6) 

_i , 

As we have Euclidean space near the axis of symmetry, we can set the radius r = R e 2 , 
with r the radius in Euclidean geometry. The Euclidean circumference turns out to be 

C = 2vrr = 2nR e"K (B.7) 

for a fixed value of R. In case that R is variable, i.e. it is no more the radius of the small 
circle around the z-axis, then the radius r is given by 

r'= l R dR'e^ ir -^\ (B.8) 



for constant tp. In this case the circumference C is different: 

rR x 

C = 2ixr' = 2tt / dtfeiF-W. (B.9) 

Since the geometry near the z-axis (R — > 0) is Euclidean, because there is no source, both 
circumferences has to be equal to each other. 

lim —7=1. (B.10) 

And therefrom 

i? e"^ 

lim fiZ_ =- = 1. (B.ll) 



i„ 1 

With the help of the method by "de l'Hopital" , we get for the limes 

l + lR^Bip(T,R) 1 
lim 2 , dR ^ K - ' = 1 ^ = 1. B.12 

e ^r(T,ii) e F(T,o) 

Obviously the last equation holds only if 

T(T,0) = 0, VTgR. (B.13) 
And especially it follows of course 

T(0,0)=0. (B.14) 

The derivative of (IB. 131) with respect to the time T is clearly zero, too. Then the assertion 
has been proved: 

T(T,R)= / dR'T, R ,(T,R'). (B.15) 
Jo 

QED 

Finally, we see that, since T,r is the energy density, the function T(T,R) describes the 
energy of the massless scalar field ip (up to a factor 2tt) contained in a disk with radius R 
and a thickness of Az = 1 at the time T. 



B.l Computation of T(T,R) 
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B.l Computation of T(T, R) 

We evaluate now the energy T(T, R) explicitly. Therefore we integrate the left side of 
(IB. 151) with respect to the radius R. 

pR pR 

r(T, R) = / dR! T, R , (T, R!) = \ / dR! R! (^, 2 T +^, 2 R/ ), (B.16) 
Jo Jo 



with 



POO 

ijj(T, R)= dkk J (kR) [a{k) e ikT + a*{k) e~ ikT ] (B.17) 
Jo 



The dependence T(T,R) on R is due to the Bessel function in tf)(T,R). Indeed only the 
Bessel function is subject to the radius R. So the integration over R concerns this function 
only. In [5] there is the formula 11.3.29 for integrating this kind of Bessel function product. 
The results are: 



I 



R D 

dR!R! UkR!) Uk'R!) = — ^ \k J^kR) J^k'R) - k' J (kR) Ji(k'R)] 



o 

R 



k 2 - k> 

[ dR'R' J x {kR!) h{k'R') = - - /2 [k h(kR) J^k'R) - k' h{kR) h(k'R)] 
Jo k k' 

(B.18) 

Therewith the integration over R has been performed, and the resulting large expression 
can be cast in a compact form, if the following is set: 

ai(k) = a(k) ai(k') = a(k') 

a 2 (k) = a*{k) a 2 {k') = a*{k'). (B.19) 

The integration over R produces eight addends. It is advantageous to notice that four of 
them are nothing but the complex conjugate of the other four. Together with the definitions 
in ( 1B.19I) . this property is helpful to cast the function T(T, R) in a very short form: 

POO POO 

T(T,R) = \\ dk dk'[K n ai(Jfe) ai(fc') +K* U a 2 (k) a 2 (k') 
Jo Jo 

+K 12 ai {k) a 2 {k') + K* 12 a 2 {k) ai{k')] , (B.20) 



where 



K* u = K 22 

KJ 2 = K 21 , (B.21) 
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with 

K n = Kn[Jfc, k'- T,R} = k k' - R 2 e - (fc+fe ' )T \k h{kR) J^k'R) - k' J^kR) J 2 (k'R) 

k 2 — k' 

-k h(kR) J {k'R) + k' J (kR) h(k'R)] , 

K12 = K 12 [k, k'- T,R} = k k' - - /2 e' {k ' k ' )T [k h{kR) J (k'R) - k! 3 (kR) J 1 {k'R) 

+k J 2 {kR) Ji{k'R) - k' h(kR) 3 2 {k'R)]. 



Appendix C 



Ricci Tensor and Ricci Scalar of the 
Einstein-Rosen Metric 

In this appendix we present the calculation of the Ricci tensor and the Ricci scalar of the 
Einstein-Rosen metric. The computation has been performed with the help of the software 
MAPLE 6, in particular with a specific package, "with(tensor)", wherein the two tools 
"RICCI" (for Ricci tensor) and "RS" (for Ricci scalar) are contained. In the following two 
sections there are the results of both computations. The motivation for this task are shown 
in the section 1.4. 

The Einstein- Rosen metric (12. 32ft . whereof the Ricci tensor and the Ricci scalar have been 
elicited reads: 

/ _ e r(T,R)-f(T,R) o \ 

e r(T,fl)-v(T,fl) o 

g*P= o o R 2 e -i> o 

\ e^ / 



C.l Ricci Tensor 



By virtue of the symmetry of the Einstein-Rosen metric the Ricci tensor has the following 
five non-vanishing components: 



R, 



1 



00 



— [R T, RR -R ip, RR -R T, TT +R i/j, tt +T, R -ij}, R -R i\)? T ] 



Roi 
Rio 



r> 

R 

r> 

R 



Rll = ~2R [ R ^ ,RR ~ R r,RR +R T,TT ~ R ^' TT + ^ ,R +T,R ~ R 1 



(C.l) 
(C.2) 

(C.3) 
(C.4) 
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R 22 = ~\R [R Tp,RR +1p,R ~R IpiTT ] (C.5) 

R33 = ^ e^~ r [R ^, RR -R ip,TT +M- (C.6) 

We remark that the two components R i and Ri are equal to each other. This is rather 
what we expected, since the Ricci tensor is symmetric in the two indices. 

C.2 Ricci Scalar 

For the Ricci scalar we obtained the following expression: 

R = ~~2R e ~ F+f t 2i? r ' TT ~ 2R T ' RR +2R ^' RR ~ 2R ^' TT +2 +R ^ ~ R ] • 

(C.7) 
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